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Abstract 

Explicit solutions for ballistic aggregation of dust-like 
matter, whose particles stick inelastically upon colli- 
sions, are considered. This system provides a model 
of large-scale structure formation in cosmology within 
the Zel'dovich approximation. In particular we show 
the equivalence of two different representations of so- 
lutions proposed in [14, 5] for a flat ID flow, extend 
these representations to cylindrically or spherically 
symmetric flows, and provide explicit counterexam- 
ples showing how exactly these representations break 
down in the case of non-symmetric flow. 



1 Introduction 

We consider here explicit solutions for ballistic aggre- 
gation of dust-like matter, with particles that stick 
upon collisions absolutely inelastically. This system 
may be considered as a model for the large struc- 
ture formation in the Universe within the Zel'dovich 
approximation [16]. First we briefly recall how this 
model comes about. 

Consider a flat Einstein-de Sitter universe contain- 
ing massive dust-like matter that moves in a self- 
consistent gravitational field. This matter partici- 
pates in two types of motion, the homogeneous cos- 
mological Hubble expansion and the evolution of lo- 
cal perturbations of density and velocity, which leads 



* Corresponding author: ansobol@obs-nice.fr 
1 Physics Department, M. V. Lomonosov Moscow State Uni- 
versity, Russia 

2 Radiophysics Department, N. N. Lobachevskn Nizhny 
Novgorod State University, Russia 

3 Observatoire de la Cote d'Azur, Nice, France. 

international Institute for Earthquake Prediction Theory 
and Mathematical Geophysics of the Russian Academy of Sci- 
ences, Moscow, Russia 



to formation of large-scale inhomogeneities. 

To describe the latter kind of dynamics in co- 
ordinates comoving with the Hubble expansion, 
Ya. B. Zel'dovich proposed [16] to use a nonlinear 
function of time, the perturbation growth factor, in- 
stead of the cosmological time. In the new time, the 
gravitational force turns out to be approximately bal- 
anced by the 'Hubble drag' caused by the continuous 
expansion of the spatial scale. Within this approxi- 
mation fluid elements move ballistically, and the evo- 
lution of local perturbations is described by the fol- 
lowing system of equations: 



| + V.(^)=0, 



d(pu) 

at 



V • (jni® u) = 0, 



(1) 
(2) 



expressing conservation of respectively mass and mo- 
mentum. Here t is the new time parameter, p(x,t) 
is the mass density field, u(x,t) is the velocity field, 
x is the spatial position comoving with the Hubble 
expansion, and u £g) u denotes the tensor with com- 
ponents UiUj. 

When trajectories of fluid elements cross, the ve- 
locity field ceases to be single- valued, and several 
streams passing through the same spatial location 
form. Numerical simulations show (see e.g. [11, 15]) 
that the gravitational interaction between streams 
confines the matter to domains of relatively small 
width and high density (Fig. 1 on page 2), causing 
it to aggregate in wall-like, filament-like, and cluster- 
like structures. This process is called ballistic aggre- 
gation. 

In [7, 8] (see also the monograph [9] and [6]) the 
Burgers equation 

du . . 
+ (u ■ V)u = pAu 
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Figure 1: Matter concentration in regions of multi- 
streamed flow (2D numerical simulations of [11]). 

was proposed to describe ballistic aggregation. In- 
deed, in the vanishing viscosity limit fi — > solu- 
tions to the Burgers equation develop singularities in 
the form of infinitcsimally thin walls, filaments, and 
point clusters, whose shape is tantalizingly similar to 
the large-scale cosmological structure. However this 
agreement is limited, as the Burgers equation is writ- 
ten in terms of the velocity field and fails to account 
for conservation of mass or momentum. Quantita- 
tively, direct iV-body simulations of self-gravitating 
matter disagree with solutions of the Burgers equa- 
tion for identical initial data. 

These limitations can be overcome in the case of 
ID flow with flat symmetry (note that in the ID case 
a thorough analytical study is possible even without 
the aggregation approximation [1]). Several authors 
[14, 10, 5, 3] independently introduced exact solu- 
tions for ballistic aggregation in a ID flow of particles 
that stick upon collisions with exact conservation of 
mass and momentum. In particular, in [14] and [5] 
these solutions were expressed in the form of vari- 
ational principles, so that the density and velocity 
fields at any given time could be constructed by min- 
imization of certain functionals. 

It is the purpose of the present work to explore 
the possibility of extending these results to the case 
where there is no flat symmetry. In Section 2 we 



prove the equivalence of the variational principles 
proposed in [ERS] and [S] 1 in the ID case with flat 
symmetry. To make the exposition self-contained 
we recall the necessary details of constructions of [S] 
and [ERS] , in particular because Ref. [S] is rather dif- 
ficult to obtain. In Section 3 the variational principles 
are extended to the case of cylindrical and spherical 
symmetry. 

Both variational principles [S] and [ERS] can be 
extended formally to the case of 2D or 3D non- 
symmetric flow, but the corresponding constructions 
no longer coincide. Moreover, there exist coun- 
terexamples showing that neither variational princi- 
ple is valid in its extended version already for non- 
symmetric 2D flows. The multidimensional exten- 
sions of the variational principles [S] and [ERS] and 
explicit counterexamples are presented in Section 4. 

The second and third authors are grateful to the 
French Ministry of education for its support. The 
work of SG was also supported by RFBR (grant 05- 
02-16517) and the Russian Leading Schools in Sci- 
ence support program (NSh-5200.2006.2). The work 
of AS was supported by RFBR (grant 05-01-00824). 

2 Ballistic aggregation in ID 

2.1 Clusters and free particles 

In the case of flat symmetry the density field p(x, t) 
and the velocity field u(x, t) both depend on a sin- 
gle coordinate. Mass concentration due to inelastic 
collisions results in the formation of parallel mas- 
sive 'walls' in three-dimensional space, which can be 
considered as point-like structures (clusters) in the 
x axis. We consider first the evolution of a purely 
discrete mass distribution. 

Let N particles in the x axis be located initially 
(i.e., at t — 0) at points y\ < y<i < • • • < yjy and have 
velocities u\, u?, . . . , it at and masses mi, mi, . . . , mpj . 
As time passes, i-th particle describes the trajectory 
Xi(t) = yi + U{t unless it collides with another parti- 
cle. After an inelastic collision there forms a cluster 
containing the whole mass of the group of particles 

1 In the sequel, [14] and [5] are referred to as [S] and [ERS] 
correspondingly. 
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that stuck together. Since trajectories of free parti- 
cles cannot cross, such clusters will always contain 
groups of particles numbered contiguously, i.e., with 
i~ < i < i + . Mass, velocity and coordinate of each 
cluster coincide with the total mass and center-of- 
mass velocity and coordinate of the corresponding 
group of particles, defined for the free motion: 



m 



= V mi, u = — V 

4+ i <i<i + 

— 12 {yi + Uit)mi. 



UiTTLi, 



-<i<i+ 



(3) 



-<i<i+ 



In order for i -th particle to stay free at time t > 0, 
its trajectory must not cross center-of-mass trajecto- 
ries of any groups of particles immediately adjacent 
to it. Therefore for any i', i" such that 1 < i' < 
io < i" < N, the coordinate Xi (t) must satisfy the 
inequalities 



m 



—jYl (f* + u ^) m ^ < ^oC*) = Via + u to~t 
m' t * 



i'<i<i 



(4) 



rn. 



where 



m 



12 mi 

i'<i<io 



m 



E 

i a <i<i' 



mi. 



ment po(y) dy. Equation (4) then assumes the form 



Jy°(y + u (y)t)po(y)dy 



f y m Po(y)dy 



<yo + u a (ya)t 



< 



JX" Po(y)dy 



(5) 



If these inequalities are satisfied for any y' < y < y" 
and at least one of them is strict, the particle located 
initially at yo will stay free at t > and will be 
located at x(y , t) = y + u (yo)t. 

On the other hand, let a group of particles located 
initially at y~ < y < y + stick together at time t into 
a cluster surrounded by free particles. Then mass, 
velocity, and coordinate of this cluster are given by 



rn 



Po(y) dy, u = — / u (y)p (y) dy, 
mj y - 



l -(\y 

m J y - 



(6) 



+ u (y)t)p (y)dy. 



For any particle specified by the initial position y, 
equations (5) and (6) together give its position x(y, t) 
at time t > without invoking the dynamics at in- 
termediate times. 

There are two ways to construct an explicit formula 
for x(y,t) from (5) and (6), which are discussed in 
the following two subsections. Both of them involve 
minimizing suitable functionals and can therefore be 
considered as (generalized) variational principles. 



Using these conditions and equations (3), one can find 
coordinates and masses of all particles and clusters 
at t > without explicitly computing dynamics of 
the system at intermediate times. Thus equations 
(4) and (3) give a complete solution for a discrete 
initial mass distribution. 

In the case of a continuous mass distribution the 
initial density and velocity fields are given by the 
functions po(x) and uq(x). In this case one replaces 
summation in equations (3) and (4) with integra- 
tion, the index i with the continuous initial coordi- 
nate y, and mass of i-th particle with the mass ele- 



2.2 The variational principle [ERS] 

Introduce in (5) and (6) a 'mass coordinate' 



m{y) 



-f 



Po(v) dri 



(7) 



and let 



$ (m) 



/y(m) 
VPo(v)dv, Uo(m) 



/y(m) 
Mv)Pa(v) dv, 



$ t (m) = $ (m) + tU (m) 



(8) 
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The lower limits of integration in (7), (8) may be 
chosen arbitrarily. It can easily be checked that the 
inverse function y(m) and Uo(y) have the form 



y(m) 



d$o(m) 



u (y(m)) 



dUo(m) 



(9) 



dm ' dm 

Hence, as long as some particle is not absorbed into 
a cluster, its coordinate is 

rf$o d(tU ) d$ t 

y+,io(j/)t = ^ + ^r = ^- (10) 

However after crossing of trajectories this represen- 
tation of x(y, t) breaks down. To give an expression 
for x(y, t) that is valid for all times, it is convenient 
to consider the cases of a free particle and a cluster 
separately. 

In the new variables, condition (5) determining if 
a particle stays free at time t > takes the following 
form: if mo = m(y ), then for any m' < m < m" 



$ t (m ) - $t(m') $ t (m") - <& t (m ) 



mo — ml 



< 



m" 



m 



(11) 



In other words, at values of m corresponding to free 
particles the function <f> t (m) coincides with its convex 
hull conv$ t , i.e., the maximal convex function not 
exceeding $ t (Fig. 2.2). One can visualize the graph 
of a convex hull as the shape of an clastic thread 
tightly wrapped around the graph of <& t from below. 
Since the function $ t is differentiable, its convex hull 
is also differentiable, and at common points of their 
graphs the derivatives d$ t /dm and cZ(conv $ t )/dm 
coincide. Therefore at these points equation (10) 
takes the form 



x(y,t) 



d$ t _ rf(conv$i) 
dm dm 



Conversely, if by the time t a group of particles 
initially located in the segment y~ < y < y + forms 
a cluster surrounded with free particles, then condi- 
tion (11) is violated at all internal points of the corre- 
sponding segment m~ = m(y~) < m < m+ = m(y + ) 
and the graph of $ t (m) lies above the chord connect- 
ing 3> t (m _ ) and <E> t (m + ), and thus above the convex 
hull (Fig. 2.2). Moreover, (6) implies that 

$ t (m+) - $ t (m~) d(conv<I> t ) 

m+ — m~ dm 



mo m 



mi 



m 



+ 



m 



Figure 2: Construction of a solution according to the 
variational principle [ERS]. At mo condition (11) is 
satisfied with any m! < m < m", while at mi it 



is violated (e.g. for m' = m 



the slope 



of the left dotted line is positive while that of the 
right one is negative). On the segment (m~,m + ), 
the convex hull conv $t(m) (dashed line) does not 
coincide with <!> t (m). Thus particles corresponding 
to points of the segment (m~, m+) stick together and 
form a cluster of the mass m + — m~ by the time t. 



Thus in clusters as well as in intervals of continuous 
mass disribution the map x(y, t) is determined by the 
derivative of the convex hull of $4: 



x { y { m),t) = <™2 tim)) 



(12) 



In [ERS], this expression is called the Generalized 
Variational Principle by analogy with the Hopf-Lax- 
Olcinik variational principle used to construct solu- 
tions to the Burgers equation. The latter principle 
itself is a variant of the least action principle in me- 
chanics. In the present paper, (12) will be referred to 
as the variational principle [ERS]. Analogous repre- 
sentations for x(y, t) were introduced independently 
in [10] and [3]. 

In the case of a constant initial density, when 
po(y) = 1 and m = y, the variational principle 
[ERS] coincides with the Hopf-Lax-Oleinik varia- 
tional principle. 
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2.3 The variational principle [S] 

Another representation for x(y,t) was introduced a 
decade earlier by A. I. Shnirel'man in [S]. To simplify 
notation, introduce the displacement field 

= x(y,t) - y. 

As long as trajectores of particles do not cross, 
£t{y) = v>o(y)t holds and the coordinate x(y,t) = 
V + €t(y) is monotonic as a function of y: 



y' + Uv')<y" + tt(y") 



(13) 



whenever y' < y" . 

We consider below only displacement fields for 
which the integral J £%{y)po{y)dy exists. This is the 
case if the initial density po(y) and velocity Uo(y) 
(and therefore the displacement field) either decrease 
at infinity fast enough or are periodical. 

Call a displacement field feasible if trajectories of 
particles do not cross, i.e., if condition (13) is satis- 
fied for any y' < y" . According to the variational 
principle [S], the solution x(y, t) is given by 

x(y,t) = y + it(y), 

where £t (y) is a feasible displacement field that mini- 
mizes the following norm of discrepancy with respect 
to uo(y)t: 



kt(y) - u (y)t\ pa(y)dy 



(14) 



In the simplest case when trajectories of particles 
do not cross, |t(y) — Uo(y)t; this displacement field is 
feasible and therefore the discrepancy vanishes. After 
crossing of trajectories a feasible displacement field 
minimizing (14) can no longer coincide with u$t. 

Mathematically, minimization problem (14) is 
equivalent to the orthogonal projection of the field 
uo(y)t to the set of feasible displacements, which we 
denoted hereafter by X. Here orthogonality is under- 
stood in the sense of the functional scalar product 



v = j ' Z(y)v(y) Po(y) dy, 



and the square of the norm of discrepancy (14) is its 
scalar poduct with itself: 

IIC* — -"o*|| 2 = J |6(y) - tu a (y)\ 2 p (y) dy 

= (6 - uot) • (6 - not). (15) 

Observe that feasible displacements form a convex 
set: if displacement fields £i(y), £2(2/) are feasible, 
then so is their 'mixture' (1 — a)£i(y) + 0^2(2/) for 
any < a < 1. Moreover, a set of feasible displace- 
ments is closed: a limit of any sequence of feasible 
displacements is also feasible. It is a well-known fact 
of functional analysis (see, e.g., [2]) that an orthogo- 
nal projection to a closed set X exists and is defined 
uniquely. Therefore the variational principle [S] al- 
ways gives a uniquely defined result. 

We now show that this result coincides with the 
result given by the variational principle [ERS]. De- 
note the displacement field defined by the map x(y, t) 
from (12) by 

- rf(conv$ t (m)) 

6(j/M) = ^ y(m). 

It suffices to check if this displacement field min- 
imizes (14), i.e., if for any other feasible displace- 
ment £ in X 

1 1 2 ~ 1 1 2 

||£ - uot 1 1 > ||& - itot|| • 

It is clear from Fig. 3 that the orthogonal projection 
to the closed set X must satisfy an even stronger 
inequality, 



1 1 2 ~ 2 

uot\\ + M - 6 , 



(16) 



and it is this condition that we will check. 
Using (15) rewrite inequality (16) in the form 

(I* -0 ■(«(>*- &) >° 

and compute its left-hand side: 

(it - ■ (u n t - It) 
= j (it(v)-m)(Mv)t-it(v))po(v)> 



(17) 



(18) 
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Figure 3: Orthogonal projection to a convex set X 
in a Hilbert space. Shown is a section of the set X 
with a hyperplane passing through the elements uot, 
£t and £. Since the point £ t is closest to u t in the 
set X, this set lies outside the sphere determined by 
the radius (u t,^ t ). Therefore in the triangle con- 
necting uot, £t, and £ the angle at £t must be obtuse, 
i.e., the square of the corresponding segment is larger 
than the sum of squares of two other segments (in- 
equality (16)). 



Observe that 



wo* -It = y+u t-x(y,t) 



<9$ t <9(conv<J> t ) 



dm 



dm 



, (19) 



where y — y(m). This expression is nonzero on 
a set of segments (m^~,m^~) of the m axis where 
3>t(m) > conv<j> t (m) is satisfied. Therefore (18) 
decomposes into a sum of integrals over segments 
(y(mi ),y(m+)) = (y~,y+), on each of which x(y,t) 
is constant, 



x(y(m),t) 



m 



and the integrand in (18) is equal to the product of 
the continuous function (uo(y)t—£t(y))po(y) with the 
function 

it(y) - Z(y) = Xi - (y + £(y)), 

which decreases because £(y) is a feasible displace- 
ment field. According to Bonnet's formula for the 
average value of the product of a monotonic and a 



continuous function 2 , there exists a point y% in the 
segment ,y^~) such that 



(xi-y- £(?/)) (u (y)t - it(y))po(y) dy 
= {x i -y i -£(^ )) _(u (y)t-tt(y))po{y)dy 

- (xi - yf - i{y+)) [ V \u (y)t - £ t (y))p (y) dy. 



(20) 



Now take into account that at y i < z <yf 

I _{u (y)t - it{y))po{y) dy 

r m(z) F)(rnn v (bA\ 

dm 



-L 



' :) (d$ t 9(conv^) 
\ dm dm 

= $ t (m(z)) - conv$ t (m(z)) > 0, (21) 



and this integral vanishes only when z = y\ . There- 
fore 



(ua(y)t - it{y))po{y)dy 

= - f y \u (y)t - £ t {v))po(v) dy, (22) 

•*Vi 



and finally 



.,7 



(xi-y- £(y)){u (y)t - £t(y))po(y) dy 
-(^ + + e(y 4 + )-^ r -C(yD)x 

x ($t(m(j/i)) - conv$ t (m(j/i))) > 0. (23) 



Thus the whole integral (18) is positive as a sum of 
positive terms. Consequently, the displacement field 
£t satisfies the inequality (17) and minimizes the vari- 
ational principle [S] (14). Conversely, the unique op- 
timal displacement field in (14) must minimize the 
variational principle [ERS]. 

2 If the function f(x) is monotonic and g{x) is continuous, 
then there exists a point c in the segment (a, ft) such that 
Jaf(x)9(x) dx = f(a) J°g(x) dx + /(ft) J*g(x) dx. 
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3 Ballistic aggregation in cylin- 
drically and spherically sym- 
metric cases 

Let the initial density p and velocity u fields de- 
pend on the (i-dimensional vector y and be symmetric 
with respect to the origin: 



po(y) = po(\y\), u (y) 



M\v\ 
\y\ 



-y- 



For d = 2 this symmetry is cylindrical and for d = 3, 
spherical. The results of this section hold for any 
d > 1, including the nonphysical case of higher di- 
mensions. 

To extend the construction of the variational prin- 
ciple [ERS] to this case one has to replace the mass 
element po(y) dy by po(\y\)fld dy in (5) and (6). Here 
ild = 27r d / 2 /T(d/2) is the (hyper)surface area of the 
unit sphere in d-dimensional space. Arguing similarly 
to Section 2.1, introducing the mass coordinate 



m(\y\) 



\v\ 



po(r))^ddr), 



and setting 



$o(m) = / r]p (r])ri d dri, 
Jo 

ry{m) 

U {m) = u a (r])po(r])n d dr), 
Jo 

$ t (m) = $ Q (m)+tU (m), 



where 



y(m) 



rfgo(m) 
dm 



u {y(m)) 



dU {m) 
dm 



we get an explicit solution that formally coincides 
with (12), but with a new definition of the func- 
tion $ t : 



4 On constructing solutions in 
the non-symmetric case 

Can the above explicit construction of solutions to 
the ballistic aggregation problem be extended to the 
case of non-symmetric flows? 

An existence theorem for a solution to the ballis- 
tic aggregation problem that respects the mass (1) 
and momentum (2) conservation laws was proved in 
[13]. However, this proof is non-constructive and thus 
gives little information about the structure of the so- 
lution and provides no method for approximating it 
numerically. 

On the contrary, both variational principles [S] and 
[ERS] provide structural information and approxi- 
mation procedures. Moreover, for a potential ini- 
tial velocity field u (y) these variational principles 
have natural multidimensional extensions requiring 
no specific assumptions on symmetry of the flow. 
However it is possible to construct examples of asym- 
metric flows for which these extensions give physi- 
cally incorrect answers, as we show below. 

4.1 A multidimensional extension of 
the variational principle [ERS] 

To extend the variational principle [ERS] (12) to the 
multidimensional case, introduce a vector "mass" co- 
ordinate m of the same dimension as the vector y, 



dm 



dy 



= po(y), 



and define functions $ ( rn ), Uo{ m ) such that ana- 
logues of equations (9) hold: 

y(m)=V m $ , u (m) = X7 m U . 

Then it is natural to define x(y,t) similarly to (12): 

x(y(m), t) = V m conv($ t + U t). (24) 

Whatever the way in which m(y), $o(m) 
and Uo(m) arc defined in general, for the uniform 
initial density po(y) = 1 it suffices to set 



m = y, $ (m) 



m 



(25) 
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and take Uq for the initial velocity field potential Uq. 
Still even in this case it is possible to show that 
the proposed modification of the variational princi- 
ple [ERS] leads to incorrect results. 

Consider a spherical wave expanding fom the ori- 
gin. Suppose that initially matter is at rest while the 
velocity potential Uq is constant everywhere except a 
small area about the origin. Asymptotically for large 
times the size of this area and details of the veloc- 
ity field inside it become inessential and the initial 
velocity potential can be approximated with 



u (y) 



o. 



y = o, 



U>0, \y\>0. 



At time t > the radius R(t) of the spherical wave is 
defined by the derivative d(conv $ t )/d\y\ at \y\ = 0, 
where 



My)= lJ Y + u (y)t, 



conv $ t (y) 




< \y\ < VWt, 
2Ui< \y\ . 



(26) 



So the proposed extension of the variational princi- 
ple [ERS] predicts that the radius of the spherical 
wave grows as V2Ut regardless of the dimension of 
space. 

On the other hand, the full mass of the wave front 
at time t is equal to the mass of matter initially dis- 
tibuted over the area contained within the wave front, 
i.e., proportional to R d , where d = 2 in the cylindri- 
cal case and d = 3 in spherical case. Since matter 
moves only in radial directions, the full momentum 
is preserved in any cone with the vertex at the origin. 
Hence the product of the velocity of the wave front 
by its mass must stay constant: R d dR/dt = const, 
so that 

R(t) cx tV(d+i). ( 27 ) 

It is this time dependence of R that appears in Sec- 
tion 3, while (26) and (24) turn out to hold only 
for d=l. 

It is interesting to determine the time dependence 
of the kinetic energy E{t). The full kinetic energy of 
the flow is the product of the mass of the vawe front 



by its squared radial velocity. Accordingly, in the 
physically correct solution (27) we have E cx R~ d cx 
t -d/(d+i) ) whilc f ormu i a (24) gi vcs E cx t( d - 2 )/ 2 ; for 

d = 3, the kinetic energy of the spherical wave turns 
out to grow indefinitely! 

Although the choice of mass coordinates m(y) is 
not uniquely determined by po alone, this freedom is 
not enough to resolve this contradiction. For a con- 
stant initial density po one can show [4] that the only 
way to introduce mass coordinates while preserving 
the spherical symmetry is (25). Thus the convex hull 
construction of the variational principle [ERS] cannot 
be extended to the non-symmetric multidimensional 
case. 

Another analysis showing incorrectness of this ex- 
tension of the variational principle [ERS] is given 
in [12]. 

4.2 A multidimensional extension of 
the variational principle [S] 

To extend the variational principle [S] to the multi- 
dimensional case, one has to define a set of feasible 
displacement fields in such a way that a solution to 
the ballistic aggregation problem would still be de- 
fined by an orthogonal projection construction. 

In one dimension a displacement field £(y) is fea- 
sible when the function x(y) = y + £(y) is mono- 
tonic, or equivalently when its primitive function is 
convex. Observe that, regardless of the dimension, 
the 'mixture' of two convex functions <&, <J/ of the 
form (1 — a)<3> + a^, where < a < 1, is always 
convex and that the differentiation operation is lin- 
ear. Therefore the class of gradients of convex func- 
tions in any dimension is a convex subset of a suit- 
able functional space. Call the displacement field 
£(y) = x(y) — y feasible if x(y) is the gradient of 
a convex function. Then a natural multidimensional 
extension of the variational principle [S] consists in 
minimizing the norm of discrepancy 

J \€t(y) - u o(y)t\ 2 Po(y) dy = min (28) 

over the class of feasible displacement fields defined 
as above. 
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Suppose that the initial vclosity field u is poten- 
tial, so that 

My) = VC/b(w), x(y,t) =y + \7U {y)t = V$ t (y), 
where 

*t(y) = ^|y| 2 + c/ (y)i. 

Evidently, at t = the function $o(y) = |y| 2 /2 is 
convex. 

As long as the Jacobian of the map y i— > x(y,t) 
is non-zero (positive), trajectories of particles do not 
cross. Since this Jacobian coincides with the determi- 
nant of the matrix of second derivatives (the Hessian 
matrix) (didj$ t ), convexity of the function <& t (y) is 
preserved. Therefore the displacement field u (y)t 
stays feasible and the discrepancy norm (28) is equal 
to zero, so the proposed form of the variational prin- 
ciple [S] gives a correct solution in this case. 

After crossing of two infinitcsimally close trajecto- 
ries the determinant of the Hessian matrix vanishes, 
the function &t(y) ceases to be convex, and the dis- 
placement field minimizing the (28) can no longer be 
given by u t. Similarly to the argument of Section 2.3 
one can show that in the spherically symmetric case 
the proposed form of the variational principle [S] (28) 
gives a correct solution. However the following exam- 
ple shows that the proposed extension gives wrong 
results in less symmetrical cases. 

Consider two identical cylindrical waves in a two- 
dimensional space with centers separated from each 
other by 2r and radii growing as R(t) = Kt 1 / 3 in 
accordance with (27). At time t > (r/K) 3 the fronts 
of these waves start to coalesce, giving rise to a flat 
portion of the wave front with matter moving longi- 
tudinally (Fig. 4). We now construct the displace- 
ment field describing the distribution of matter in 
the cylindrical portions and in the flat portion of the 
wave front for t > (r/K) 3 . 

Let the coordinates be chosen so that centers of 
waves are at (±r, 0) . At time t the flat portion is the 
segment of length 2l(t) = 2^J R 2 (t) — r 2 in the y axis. 
Matter coming to this segment from the rounded por- 
tions of the wave front sticks together and moves out- 
ward with velocity v = R ■ l/R. Observe that the 



v 




1 Kt) 


\h(*) \ 


l — r 


r J x 







Figure 4: Two cylindrical waves: an example of a 
flow for which the proposed multidimensional exten- 
sion of the variational principle [S] (28) gives an in- 
correct result. The front is shown at time t satisfying 
(r/K) 3 <t<t* = (b/2f/ 2 (r/K) 3 . 



length of the segment increases faster, at the rate 

• RR ■ R ■ I 

I = , = R ■ — > R ■ — = v. 

VR 2 - r 2 I R 

thus endpoints of the flat segment do not affect the 
motion of matter inside it. 

Let T(l) be the inverse funcion for l(t), 

T(l)^^(f + r 2 f/ 2 1 

and take I as a parameter. Since R(t) — Kt 1 / 3 , we 
have R/R = 1/3*, so that v(l) = IR/R = l/3T(l). 
Therefore at time t > T(l) the particle parametrized 
with I will be located at 

y(l,t) = l + v(l)(t-T(l)) = h + ^ ry 

The distribution of mass inside the flat segment stays 
continuous as long as the map I y(l,t) is mono- 
tonic, i.e., while dy/dl = 2/3 + (t/3)(l/T(l))' > for 
all I. This inequality is satisfied for (r/K) 3 < t < t* , 
where 

r = ^("(TOy) = G) / © 3 - 
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Here minimization is performed for I > r/y/2 because 
for < I < r/V2 the condition dy/dl > holds for 
any t. 




Figure 5: Continuous black lines starting at (±r, 0) 
show initial positions of particles located at the points 
of the wave front in the flow of Fig. 4 (for times 
(r/K) 3 <t<t* = {b/2f/ 2 {r/Kf). Note that parti- 
cles forming the flat portion of the front are initially 
located on two-segment broken lines in the central 
area of the figure. 

Observe now that initial positions of particles that 
are located at the point (0, y(l, t)) in the flat segment 
of the front at time T(l) < t < t* form a two-segment 
broken line connecting the point (0, I) with centers of 
both waves (Fig. 5). On the other hand, the preimage 
of any point in a mapping defined by the gradient of 
a convex function must itself be a convex set (a point, 
a segment, or a convex domain). Since preimages of 
some points in the displacement field are non-convex 
broken lines, this displacement field cannot be real- 
ized by the gradient of a convex function, i.e., it is not 
'feasible' in the above sense and the proposed natu- 
ral extension of the variational principle [S] provides 
an incorrect description of describes the motion of 
matter in the case of two cylindrical waves. 



References 

[1] E. Aurcll, D. Fanelli, S. N. Gurbatov, and A. Y. 
Moshkov. The inner structure of Zeldovich pan- 



cakes. Physica D, 186 (2003) 171-184. 

[2] A. V. Balakrishnan. Introduction to optimiza- 
tion theory in a Hilbert space. Lecture notes in 
operations research and mathematical systems. 
Springer- Verlag, 1971. 

[3] Y. Brenier and E. Grenier. Sticky particles and 
scalar conservation laws. SI AM J. Numer. Anal., 
35:6 (1998) 2317-2328. 

[4] L. Caffarclli and Y. Li. An extension to a theo- 
rem of Jorgens, Calabi, and Pogorelov. Comm. 
Pure Appl. Math., 56 (2003) 549-583. 

[5] E W., Y. Rykov, and Y. Sinai. Generalized vari- 
ational principles, global weak solutions and be- 
havior with random initial data for systems of 
conservation laws arising in adhesion particle dy- 
namics. Comm. Math. Phys., 177:2 (1996) 349- 
380. 

[6] S. N. Gurbatov and A. Y. Moshkov. Genera- 
tion of large-scale coherent structures in the kpz 
equation and multidimensional burgers turbu- 
lence. J. Experimental and Theoretical Phys., 
976 (2003) 1186-1200. 

[7] S. N. Gurbatov and A. I. Saichev. Probability 
distributions and spectra of potential hydrody- 
namic turbulence. Izv. Vyssh. Uchebn. Zaved. 
Ser. Radiofizika, 27 (1984) 456-468. 

[8] S. N. Gurbatov, A. I. Saichev, and S. F. Shan- 
darin. The large-scale structure of the universe 
in the frame of the model equation of non-linear 
diffusion. Mon. Not. R. Astron. Soc, 236 (1989) 
385-402. 

[9] S. N. Gurbatov, A. N. Malakhov, and A. I. 
Saichev. Nonlinear Random Waves and Turbu- 
lence in Nondispersive Media. Nonlinear Science: 
Theory & Applications. Wiley, 1992. 

[10] P. A. Martin and J. Piasecki. One dimensional 
ballistic aggregation: Rigorous long-time esti- 
mates. J. Stat. Physics, 76 (1994) 447. 



10 



[11] A. L. Melott and S. F. Shandarin. Gravitational 
instability with high resolution. Astrophys. J., 
343 (1989) 26-30. 

[12] Y. G. Rykov. On the nonhamiltonian character 
of shocks in 2-D pressureless gas. Boll. Unione 
Mat. ltd. Sez. B Artie. Ric. Mat. (8), 5:1 (2002) 
55-78. 

[13] M. Sever. An existence theorem in the large for 
zero-pressure gas dynamics. Differential Integral 
Equations, 14:9 (2001) 1077-1092. 

[14] A. I. Shnirel'man. On the principle of the short- 
est way in the dynamics of systems with con- 
straints. In Global analysis — studies and appli- 
cations, II. Lecture Notes in Math., vol. 1214. 
Springer- Verlag, 1986, 117-130. 

[15] D. H. Weinberg and J. E. Gunn. Largescale 
Structure and the Adhesion Approximation. 
Mon. Not. R. Astron. Soc. 247 (1990) 260-286. 

[16] Y. B. Zel'dovich. Gravitational instability: an 
approximate theory for large density perturba- 
tions. Astron, & Astrophys., 5 (170) 84-89. 



11 



